

I0 = 20;p = 100;
i = -100:100;r = 2*i;
j = I0./(2*pi*r.^2);
U = I0*p./(2*pi*r);
E = j*p;
hold on 
plot(j);   plot(U); plot(E); 
hold off